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ABSTRACT 

In  this  study,  the  U.S.  Navy’s  Generalized  Digital  Environmental  Model  (GDEM)  climatological  temperature 
and  salinity  data  on  a  0.5°  X  0.5°  grid  is  used  to  investigate  the  seasonal  variabilities  of  the  Japan/East  Sea 
(JES)  thermohaline  structure  and  circulations.  The  GDEM  for  the  JES  was  built  up  on  historical  (1930-97) 
136  509  temperature  and  52  572  salinity  profiles.  A  three-dimensional  estimate  of  the  absolute  geostrophic 
velocity  field  was  obtained  from  the  GDEM  temperature  and  salinity  fields  using  the  P-vector  method.  The 
seasonal  variabilities  of  the  thermohaline  structure  and  the  inverted  currents  such  as  the  Subpolar  Front,  the 
salinity  minimum  and  maximum  in  the  Japan  Sea  Intermediate  Water,  and  the  Tsushima  Warm  Current  and  its 
bifurcation  are  identified. 


1.  Introduction 

The  Japan  Sea,  known  as  the  East  Sea  in  Korea,  has 
steep  bottom  topography  (Fig.  1)  that  makes  it  a  unique 
semienclosed  ocean  basin  overlaid  by  a  pronounced 
monsoon  surface  wind.  The  Japan/East  Sea  (hereafter 
JES)  covers  an  area  of  106  km2,  has  a  maximum  depth 
in  excess  of  3700  m,  and  is  isolated  from  open  oceans 
except  for  small  (narrow  and  shallow)  straits  that  con¬ 
nect  the  JES  to  the  Pacific  Ocean.  The  JES  contains 
three  major  basins  called  the  Japan  Basin  (JB),  Ulleng/ 
Tsushima  Basin  (UTB),  and  Yamato  Basin  (YB),  and  a 
high  central  seamount  called  the  Yamato  Rise.  The  JES 
is  of  great  scientific  interest  as  a  miniature  prototype 
ocean.  Its  basinwide  circulation  pattern,  boundary  cur¬ 
rents,  Subpolar  Front  (SPF),  mesoscale  eddy  activity, 
and  deep-water  formation  are  similar  to  those  in  a  large 
ocean. 

The  JES  experiences  two  monsoons,  winter  and  sum¬ 
mer,  every  year.  During  the  winter  monsoon  season,  a 
cold  northwest  wind  blows  over  the  JES  as  a  result  of 
the  Siberian  high  pressure  system  located  over  the  East 
Asian  continent.  Radiative  cooling  and  persistent  cold 
air  advection  maintain  cold  air  over  the  JES.  The  north- 
west-southeast  oriented  jet  stream  is  positioned  above 
the  JES.  Such  a  typical  winter  monsoon  pattern  lasts 
nearly  six  months  (Nov-Apr).  During  the  summer  mon¬ 
soon,  a  warm  and  weaker  southeast  wind  blows  over 
the  JES.  Such  a  typical  summer  monsoon  pattern  lasts 
nearly  4  months  (mid-May  to  mid-Sep). 
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The  JES  thermohaline  structure  and  general  circula¬ 
tion  have  been  investigated  for  several  decades.  The 
Tsushima  Warm  Current  (TWC),  dominating  the  surface 
layer,  flows  in  from  Tsushima  Strait  and  carries  warm 
water  from  the  south  up  to  40°N  where  a  polar  front 
forms  (Seung  and  Yoon  1995).  Most  of  the  nearly  ho¬ 
mogeneous  water  in  the  deep  part  of  the  basin  is  called 
the  Japan  Sea  Proper  Water  (Moriyasu  1972)  and  is  of 
low  temperature  and  low  salinity.  Above  the  Proper  Wa¬ 
ter,  warm  and  saline  water  flows  in  through  Tsushima 
Strait,  transports  northeastward,  and  flows  out  through 
the  Tsugaru  and  Soya  Strait. 

The  TWC  separates  north  of  35°N  into  western  and 
eastern  channels  (Uda  1934;  Kawabe  1982a,b;  Hase  et 
al.  1999;  Senjyu  1999).  The  flow  through  the  western 
channel  closely  follows  the  Korean  coast  [called  the  East 
Korean  Warm  Current  (EKWC)]  until  it  bifurcates  into 
two  branches  near  37°N.  The  eastern  branch  follows  the 
SPF  to  the  western  coast  of  Hokkaido  Island,  and  the 
western  branch  moves  northward  and  forms  a  cyclonic 
eddy  at  the  eastern  Korean  Bay.  The  flow  through  the 
eastern  channel  follows  the  Japanese  coast,  called  the 
“Nearshore  Branch”  by  Yoon  (1982a-c).  More  accu¬ 
rately,  we  may  call  it  the  Japan  Nearshore  Branch  ( JNB). 
The  JNB  is  usually  weaker  than  the  EKWC.  The  TWC 
at  both  channels  decreases  with  depth.  The  EKWC 
meets  the  southward  coastal  current,  the  North  Korean 
Cold  Current  (NKCC),  at  about  38°N  with  some  sea¬ 
sonal  meridional  migration.  After  separation  from  the 
coast,  the  EKWC  and  the  NKCC  converge  and  form  a 
strong  front  (i.e.,  SPF)  that  stretches  in  a  west-east  di¬ 
rection  across  the  basin.  The  NKCC  makes  a  cyclonic 
recirculation  gyre  in  the  north  but  most  of  the  EKWC 
flows  out  through  the  outlets  (Uda  1934).  The  formation 
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Fig.  1.  Geography  and  isobaths  showing  the  bottom  topography  (in  meters)  of  the  Japan/East 

Sea  (JES). 
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of  the  NKCC  and  the  separation  of  the  EKWC  are  due 
to  a  local  forcing  by  wind  and  buoyancy  flux  (Seung 
1994).  Large  meanders  develop  along  the  front  and  are 
associated  with  warm  and  cold  eddies. 

Between  the  TWC  water  and  the  Japan  Sea  Proper 
Water,  a  vertical  salinity  minimum  (SMIN)  south  of  the 
SPF,  usually  accompanied  by  a  dissolved  oxygen  max¬ 
imum,  was  first  found  and  named  Japan  Sea  Interme¬ 
diate  Water  (JIW)  by  Miyazaki  (1952,  1953),  and  further 
depicted  by  Kajiura  et  al.  (1958)  and  Moriyasu  (1972). 
The  collocation  of  the  SMIN  and  the  dissolved  oxygen 
maximum  implies  that  the  Intermediate  Water  originates 
from  the  descending  surface  water  around  the  SPF  (Mi¬ 
yazaki  1952,  1953;  Miyazaki  and  Abe  1960).  Later,  Kim 
and  Chung  (1984)  found  a  very  similar  property  in  the 
UTB  and  proposed  the  East  Sea  Intermediate  Water  (i.e.. 


JIW  used  by  Japanese  oceanographers);  and  Kim  and 
Kim  (1999)  found  the  high  salinity  water  with  high 
oxygen  in  the  eastern  JB  (i.e.,  north  of  SPF)  and  named 
the  water  “High  Salinity  Intermediate  Water”  (HSIW; 
>34.07  psu). 

The  seasonal  variabilities  of  both  the  JES  circulation 
and  thermohaline  structure  have  been  studied  based  on 
limited  datasets  such  as  the  seasonal  sea  surface  tem¬ 
perature  (SST)  variability  (Isoda  and  Saitoh  1988,  1993; 
Kano  1980;  Maizuru  Mar.  Observ.  1997).  Based  on  sat¬ 
ellite  infrared  images  in  the  western  part  of  the  JES  and 
the  routine  hydrographic  survey  by  the  Korea  Fisheries 
Research  and  Development  Agency  in  1987,  Isoda  and 
Saitoh  (1993)  found  SST  patterns  in  winter  and  spring 
that  characterized  as  follows:  a  small  meander  of  a  ther¬ 
mal  front  first  originated  from  Tsushima  Strait  near  the 
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Korean  coast  and  gradually  grew  into  an  isolated  warm 
eddy  with  a  horizontal  scale  of  100  km.  The  warm  eddy 
intruded  slowly  northward  from  spring  to  summer.  Chu 
et  al.  (1998a)  reported  the  seasonal  occurrence  of  JES 
eddies  from  the  composite  analysis  on  the  U.S.  National 
Centers  for  Environmental  Prediction  (NCEP)  monthly 
SST  fields  (1981-94).  For  example,  they  identified  a 
warm  center  appearing  in  later  spring  in  the  East  Korean 
Bay. 

What  are  the  seasonal  variabilities  of  the  three-di¬ 
mensional  JES  circulation  and  thermohaline  structure? 
We  use  the  U.S.  Navy’s  unclassified  Generalized  Digital 
Environmental  Model  (GDEM)  temperature  and  salinity 
data  on  a  0.5°  X  0.5°  grid  to  investigate  the  problem  in 
this  study.  The  outline  of  this  paper  is  as  follows:  A 
description  of  the  GDEM  data  is  given  in  section  2.  A 
depiction  of  the  inversion  of  the  absolute  geostrophic 
velocity  is  given  in  section  3.  Explanations  of  monthly 
mean  and  anomaly  fields  are  given  in  section  4.  The 
seasonal  variabilities  of  the  three-dimensional  thermo¬ 
haline  structure  and  the  inverted  velocity  held  are  dis¬ 
cussed  in  sections  5-7.  In  section  8  we  present  our  con¬ 
clusions. 

2.  The  U.S.  Navy’s  GDEM  dataset 

The  U.S.  Navy’s  (GDEM)  global  climatological 
monthly  mean  temperature  and  salinity  dataset  is  a  four¬ 
dimensional  (latitude,  longitude,  depth,  and  month)  dis¬ 
play.  Data  for  building  the  present  version  of  the  GDEM 
climatology  for  the  JES  were  obtained  from  the  Navy’s 
Master  Oceanographic  Observational  Data  Set 
(MOODS),  which  has  136  509  temperature  and  52  572 
salinity  profiles  during  1930-97.  The  main  limitation 
of  the  MOODS  data  is  its  irregular  distribution  in  time 
and  space.  Certain  periods  and  areas  are  oversampled, 
while  others  lack  enough  observations  to  gain  any  mean¬ 
ingful  insights.  Vertical  resolution  and  data  quality  are 
also  highly  variable  depending  much  on  instrument  type 
and  sampling  expertise  (Chu  et  al.  1997a-d).  The 
monthly  distributions  of  the  JES  temperature  (Fig.  2a) 
and  salinity  (Fig.  2b)  stations  show  that  the  number  of 
temperature  stations  is  two  to  three  times  more  than  the 
number  of  salinity  stations.  January  has  the  least  profiles 
and  August  the  most.  Yearly  temperature  (Fig.  3a)  and 
salinity  (Fig.  3b)  profile  numbers  show  temporally  un¬ 
even  distribution  with  almost  no  observations  in  the 
whole  JES  in  certain  years  (e.g.,  1944,  1989  for  tem¬ 
perature,  and  1944,  1987-93  for  salinity)  and  many  ob¬ 
servations  in  other  years  (e.g.,  nearly  6500  temperature 
profiles  in  1969,  and  3700  salinity  profiles  in  1967). 
Spatial  and  temporal  irregularities  along  with  the  lack 
of  data  in  certain  regions  must  be  carefully  weighted  in 
order  to  avoid  statistically  induced  variability. 

The  basic  design  concept  of  the  GDEM  is  the  deter¬ 
mination  of  a  set  of  analytical  curves  that  represent  the 
mean  vertical  distributions  of  temperature  and  salinity 
for  grid  cells  (0.5°  X  0.5°)  through  the  averaging  of 


coefficients  for  the  curves  found  for  individual  profiles 
(Teague  et  al.  1990).  Different  families  of  representative 
curves  have  been  chosen  for  shallow,  mid,  and  deep 
depth  ranges,  with  each  chosen  so  that  the  number  of 
parameters  required  to  yield  a  smooth,  mean  profile  over 
the  range  was  minimized.  As  mentioned  by  Teague  et 
al.  (1990),  large-scale  oceanographic  features  are  gen¬ 
erally  found  to  be  similarly  represented  in  both  the 
GDEM  and  the  1982  NOAA  Climatological  Atlas  of 
the  World  Ocean  temperature  and  salinity.  The  GDEM 
appears  to  render  better  representations  of  seasonal  var¬ 
iability  and  regions  of  high  current  shear  because  of  a 
different  smoothing  method  and  a  finer  grid  spacing. 
The  GDEM  data  contains  monthly  mean  temperature 
and  salinity  (7)  S)  and  annual  mean  temperature  and 
salinity  ( T ,  S)  fields.  Interested  readers  are  referred  to 
Teague  et  al.  (1990)  for  more  information. 

3.  The  P- vector  inverse  method 

Recently,  we  developed  the  P-vector  inverse  method 
(Chu  1995;  Chu  et  al.  1998b, c;  Chu  2000;  Chu  and  Li 
2000)  to  calculate  the  absolute  geostrophic  velocity. 
This  method  contains  two  steps:  (i)  determination  of  the 
velocity  direction  and  (ii)  determination  of  the  velocity 
magnitude.  Two  necessary  conditions  for  the  inversion 
are  easily  implemented  into  this  method:  (i)  the  iso- 
pycnal  surface  does  not  parallel  the  potential  vorticity 
surface  and  (ii)  the  velocity  has  a  vertical  spiral  (Chu 
et  al.  1998b,c). 


a.  Reduced  physics 

As  pointed  out  by  Wunsch  and  Grant  (1982),  in  de¬ 
termining  large-scale  circulation  from  hydrographic 
data  we  can  be  reasonably  confident  of  the  assumptions 
of  geostrophic  balance,  mass  conservation,  adiabatic, 
and  no  major  cross-isopycnal  mixing  (except  for  water 
masses  in  contact  with  the  atmosphere).  Under  these 
conditions,  the  density  of  each  fluid  element  would  be 
conserved,  which  mathematically  is  given  by 

V  •  Vp  =  0,  (1) 

where  p  is  the  potential  density  and  V  =  in,  v,  w )  is 
the  geostrophic  velocity.  The  conservation  of  potential 
vorticity  equation  can  be  obtained  by  differentiating  (1) 
with  respect  to  z,  using  geostrophic  and  hydrostatic  bal¬ 
ances,  and  including  the  latitudinal  variation  of  the  Cor¬ 
iolis  parameter  to  give 

V  •  V<?  =  0,  (2) 

where  q  =  fdp/dz .  Equations  (1)  and  (2)  imply  that  the 
velocity  V  is  parallel  to  Vq  X  Vp. 


b.  The  P-vector  method 


Consider  the  unit  vector  P  (Chu  1995)  defined  by 


Vp  X  V<7 
IVp  X  Vgl 


(3) 
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The  existence  of  this  unit  vector  implies  the  satisfaction 
of  necessary  condition  (1). 

The  velocity,  Y  =  (u,  v,  w ),  parallels  the  unit  vector  P, 

V  =  r(x,  y,  z )  P,  (4) 

where  r  is  the  proportionality.  Applying  the  thermal 
wind  relation  at  two  different  depths  zk  and  zm,  a  set  of 
algebraic  equations  for  determining  the  parameter  r  are 
obtained 

r(tTf)  -  =  A ukm 

rmP[k)  -  r(m)P<m)  =  Avkm,  (5) 

which  are  two  linear  algebraic  equations  for  r'h  and  r(m) 
[r(i>  =  r(x,  y,  z,)].  Here 

■£)*'•  <6) 

where  p  is  the  in  situ  water  density,  and  p0  is  the  char¬ 

acteristic  value  of  the  density. 

The  determinant  of  the  coefficient  matrix  of  (5)  is  the 
sine  of  the  vertical  turning  angle  between  Pj^  and  Pj'”' 
(Chu  et  al.  1998b,c;  Chu  and  Li  2000).  The  existence 
of  solutions  of  (5)  implies  the  satisfaction  of  necessary 
condition  (2). 

For  water  columns  satisfying  the  two  necessary  con¬ 
ditions,  we  may  solve  (5)  to  obtain  r(k}  for  the  level  zk. 
There  are  N  —  1  sets  (m  =  1,2 ,  k  —  1,  k  +  1,  .  .  . , 
N)  of  equations  (5)  for  calculating  rik> .  Here  N  is  the 
total  number  of  vertical  levels  of  the  water  column.  All 
the  N  —  1  sets  of  equations  are  compatible  under  the 
thermal  wind  constraint  and  should  provide  the  same 
solution.  However,  because  of  errors  in  measurements 
(instrumentation  errors)  and  computations  (truncation 
errors),  the  parameters  r®  may  vary  with  m.  A  least 
squares  error  algorithm  is  used  to  minimize  the  error. 
For  further  details  and  validation  of  the  algorithm  pre¬ 
sented  here  see  Chu  et  al.  (1998b).  (Interested  readers 
can  obtain  the  software  by  contacting  C.  Fan  at 
fan@nps.navy.mil  or  by  visiting  our  Web  site  at  http:/ 
/www.  oc .  nps .  navy,  mil/— chu. ) 

4.  Monthly  mean  and  anomaly  fields 

The  GDEM’s  monthly  mean  data  and  inverted  ve¬ 
locity  field  show  an  evident  seasonal  variation.  The  av¬ 
erage  of  the  12  monthly  mean  fields  (7)  S)  leads  to  the 
annual  mean  field  ( T ,  S).  The  GDEM  dataset  also  pro¬ 
vides  the  annual  mean  fields.  The  monthly  minus  annual 
mean  field  is  defined  as  the  monthly  anomaly  field, 

T,n  =  T  -  T,  S,n  =  S  -  S, 

which  represents  the  seasonal  variation  relative  to  the 
annual  mean  field. 

The  seasonal  7)  S  variation  reduces  as  depth  increases. 
At  300-m  depth  there  is  almost  no  seasonal  variability. 
Thus,  we  present  only  the  horizontal  fields  at  the  surface 
and  the  intermediate  level  (150  m)  for  illustration.  In 


order  to  see  the  seasonal  variation  of  the  vertical  ther¬ 
mohaline  structure,  we  present  the  latitudinal  and  zonal 
cross  sections  of  the  monthly  mean  fields. 

To  identify  the  variability  in  both  monthly  mean  and 
monthly  mean  anomaly  fields,  we  use  a  cold  (warm) 
center  to  represent  the  temperature  minimum  (maxi¬ 
mum),  and  a  fresh  (salty)  center  to  represent  the  salinity 
minimum  (maximum). 

5.  Temperature 

a.  Sea  surface 

Although  the  monthly  SST  field  (Fig.  4)  shows  an 
evident  seasonal  variation,  the  SPF  exists  continuously 
throughout  the  year.  Its  position  is  quite  stationary,  but 
its  intensity  strengthens  in  winter  (especially  the  east 
part)  and  weakens  in  summer.  Such  a  pattern  is  similar 
to  an  earlier  description  (Maizuru  Mar.  Observ.  1997). 
The  location  of  the  SPF  in  spring  is  quite  consistent 
with  Isoda  and  Saitoh’s  (1988)  estimations  using  ten 
NOAA-8  satellite  Advanced  Very  High  Resolution  Ra¬ 
diometer  images  in  spring  1984. 

The  SST  is  found  always  higher  than  5°C  year-round 
in  the  UTB  and  the  YB,  consistent  with  Kim  et  al.’s 
(1999)  observational  studies.  The  maximum  SST  gra¬ 
dient  is  found  as  16°C/100  km  near  137°E,  40°N  in  Feb¬ 
ruary  and  March,  and  the  minimum  SST  gradient  is 
found  as  8°C/100  km  from  July  to  September.  The  SST 
gradient  across  the  SPF  is  two  times  as  strong  in  winter 
as  in  summer.  The  weakening  of  the  SPF  in  summer  is 
caused  by  the  faster  warming  of  the  water  mass  north 
of  rather  than  south  of  the  SPF  in  spring.  North  of  the 
SPF  a  second  front  occurs  (bifrontal  structure)  during 
the  fall-to-winter  transition  season,  especially  in  No¬ 
vember  and  December.  This  front  parallels  the  Russian 
coast  with  the  maximum  SST  gradient  around  4°C/100 
km  in  November. 

The  monthly  mean  SST  anomaly  field  obtained  from 
the  GDEM  data  (Fig.  5)  shows  a  basin-scale  cold  anom¬ 
aly  from  December  to  May  (winter  monsoon)  and  a 
basin-scale  warm  anomaly  from  June  to  October  (sum¬ 
mer  monsoon).  The  strongest  basin-scale  cooling  occurs 
in  February  with  the  coldest  anomaly  center  (7’.ln  < 
— 9°C)  appearing  in  the  central  JB  (40°-42°N,  135°- 
137.5°E);  and  the  strongest  basin-scale  warming  occurs 
in  August  with  the  warmest  anomaly  center  (7’m  >  1 1°C) 
appearing  near  Peter  the  Great  Bay  (PGB). 

The  typical  winter  anomaly  pattern  is  featured  by  (i) 
all  negative  values  and  (ii)  eastward  expansion  of  a  cold 
anomaly  from  the  PGB.  After  the  winter  monsoon  onset, 
the  cold  SST  anomaly  first  occurs  in  the  north  to  north¬ 
western  boundary  from  Tartar  Strait  to  PGB  in  Novem¬ 
ber  with  a  minimum  value  around  —  1°C.  In  December, 
the  cold  anomaly  sweeps  the  whole  JES  basin  with  the 
coldest  center  (— 4°C)  at  the  north  to  northwestern 
boundary.  The  PGB  cold  center  becomes  evident  with 
a  minimum  value  of  — 6°C  in  January.  The  PGB  cold 
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Fig.  2.  Spatial  distributions  of  MOODS  stations  during  1930-97:  (a)  temperature  and  (b)  salinity. 
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Fig.  3.  Temporal  distributions  of  MOODS  stations  during  1930- 
97:  (a)  temperature  and  (b)  salinity. 


the  year.  South  of  the  SPF,  the  temperature  changes  from 
5°C  to  9°C.  The  SPF  meandering  at  131°E,  134°E,  and 
138°E  forms  several  mesoscale  eddies.  The  SPF  me¬ 
andering  near  Okin  Gunto  (134°E)  in  spring  was  pre¬ 
viously  reported  by  Isoda  and  Saitoh  (1988,  1993).  Be¬ 
sides,  a  second  front  occurs  (bifrontal  structure)  south 
of  the  SPF  along  the  west  coast  of  Japan  during  the 
winter  and  spring  seasons. 

Monthly  mean  temperature  anomaly  Tm  (Fig.  7) 
shows  (i)  a  north-south  asymmetric  pattern  with  weaker 
seasonal  variability  north  of  rather  than  south  of  the 
SPF,  and  (ii)  weakening  of  the  SPF  in  the  winter  and 
strengthening  of  the  SPF  in  the  summer. 

An  evident  cooling  occurs  south  of  the  SPF  from 
January  to  May  (winter  monsoon)  with  a  large  cold 
anomaly  center  (Tan  <  —  2°C)  appearing  in  the  YB  (36°- 
39°N,  132.5°-136°E)  from  March  to  April.  Much  weak¬ 
er  cooling  appears  north  of  the  SPF  during  the  same 
period  (Jan-May).  Such  a  differential  cooling  causes 
the  northward  increase  of  7’m  across  the  SPF,  which 
implies  the  reduction  of  the  SPF  strength. 

An  evident  warming  occurs  south  of  the  SPF  from 
June  to  October  (summer  monsoon)  with  a  large  warm 
anomaly  center  (T.m  >  2°C)  appearing  in  the  YB  from 
August  to  September.  Much  weaker  warming  appears 
north  of  the  SPF  during  the  same  period  (Jun-Oct).  Such 
a  differential  warming  causes  the  northward  decrease  of 
7’.m  across  the  SPF,  which  implies  the  enhancement  of 
the  SPF  strength. 


center  expands  eastward  to  138°E  and  occupies  the  mid- 
JES  north  of  the  SPF  (39°-43°N)  with  the  coldest  center 
of  —  9°C  in  February.  The  spring  (Mar-May)  pattern  is 
featured  by  weakening  of  the  mid-JES  cold  center  A 
typical  summer  anomaly  pattern  is  featured  by  (i)  all 
positive  values  and  (ii)  eastward  expansion  of  a  warm 
center  from  PGB.  After  the  summer  monsoon  onset,  a 
strong  SST  warm  center  (8°C)  occurs  near  PGB  in  July. 
The  PGB  warm  center  strengthens  and  expands  eastward 
to  138°E  and  occupies  the  mid-JES  north  of  the  SPF 
(39°-43°N)  with  the  warmest  center  of  11°C  in  August. 
The  autumn  (Sep-Oct)  pattern  is  featured  by  weakening 
of  the  mid-JES  warm  center.  The  fact  that  the  coldest 
(warmest)  center  occupies  the  mid-JES  north  of  the  SPF 
(39°-43°N)  in  winter  (summer)  suggests  that  the  SST 
seasonal  variability  is  greater  in  the  north  than  in  the 
south  of  the  SPF.  This  result  is  consistent  with  an  earlier 
analysis  on  the  NCEP  SST  data  (Chu  et  al.  1998b). 

b.  Intermediate  level  (150  m) 

The  seasonal  thermal  variability  at  the  150-m  depth 
is  much  weaker  than  that  at  surface  (Fig.  6).  The  SPF 
still  exists  throughout  the  year  and  is  located  at  almost 
the  same  location  as  at  the  surface.  North  of  the  SPF, 
the  temperature  is  uniformly  cold  (1°-3°C)  throughout 


c.  Zonal  cross  sections  (37°  and  43°N) 

The  zonal  cross  sections  (37°  and  43°N)  of  monthly 
mean  temperature  show  a  strong  seasonal/permanent 
thermocline  structure  south  of  the  SPF  (Fig.  8a)  and  a 
strong  seasonal/weak  permanent  thermocline  structure 
north  of  the  SPF  (Fig.  8b). 

South  of  the  SPF  at  37°N  (Fig.  8a),  a  permanent  ther¬ 
mocline  is  located  at  80-125  m  and  appears  year-round 
with  the  maximum  strength  (0.12°C  m-1)  in  August. 
Above  it,  the  seasonal  thermocline  occurs  from  the  sur¬ 
face  to  50-m  depth  in  June  (0.15°C  m-1),  intensifies 
during  the  summer  monsoon  season  to  a  maximum  value 
of  around  0.36°C  m-1  in  August,  and  weakens  in  Sep¬ 
tember.  In  October,  the  seasonal  thermocline  erodes  and 
an  ocean  mixed  layer  (OMF)  is  formed.  In  November, 
the  OMF  is  well  established  with  the  temperature  near 
18°C  and  the  depth  around  75  m.  During  the  prevailing 
winter  monsoon  season  (Dec-Mar),  the  OMF  deepens 
to  80-130  m  with  a  westward  uplift  of  the  OMF  depth: 
80  m  near  the  Korean  coast  and  130  m  near  the  Japanese 
coast.  The  OMF  temperature  is  around  10°C.  The  OMF 
starts  to  warm  at  a  rate  of  2°C  month-1  from  March  to 
May,  and  its  depth  shoals,  respectively.  For  example, 
the  OMF  temperature  increases  from  12°C  in  April  to 
14°C  in  May,  and  the  OMF  depth  decreases  from  50- 
70  m  in  April  to  less  than  50  m  in  May.  This  process 
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Fig.  4.  Monthly  mean  temperature  (°C)  field  at  the  ocean  surface. 
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Fig.  5.  Monthly  mean  temperature  anomaly  (°C)  field  at  the  ocean  surface.  Negative  values  are  shaded. 
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Fig.  6.  Monthly  mean  temperature  (°C)  field  at  150-m  depth. 
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Fig.  7.  Monthly  mean  temperature  anomaly  (°C)  field  at  150-m  depth.  Negative  values  are  shaded. 
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(OML  warming  and  shoaling)  continues  during  the  sum¬ 
mer  monsoon  season  (Jun-Aug). 

North  of  the  SPF  at  43°N  (Fig.  8b),  the  permanent 
thermocline  is  quite  weak  and  located  near  the  surface 
to  300-m  depth  with  an  evident  westward  uplift  from 
the  Japanese  coast  (thickness  around  300  m)  to  the  Rus¬ 
sian  coast  (thickness  around  100  m).  The  seasonal  ther¬ 
mocline  occurs  from  the  surface  to  50-m  depth  in  May 
(~0.08°C  m_1),  intensifies  during  the  summer  monsoon 
season  to  a  maximum  value  of  around  0.5°C  m_1  in 
August  and  September,  and  weakens  in  October.  In  No¬ 
vember,  the  seasonal  thermocline  erodes  and  becomes 
part  of  the  permanent  thermocline,  which  weakens  dur¬ 
ing  the  prevailing  winter  monsoon  season.  In  February, 
the  permanent  thermocline  is  so  weak  that  the  water 
column  is  almost  uniformly  cold  (1°C)  west  of  136°E 
and  weakly  stratified  (<0.01°C  m_I)  east  of  136°E. 

d.  Latitudinal  cross  section  (135°E) 

The  strong  north-south  thermal  asymmetry  across  the 
SPF  is  also  identified  from  a  latitudinal  cross  section 
(135°E)  of  the  monthly  mean  temperature  (Fig.  9). 

During  the  prevailing  winter  monsoon  season  (Dec- 
Mar),  the  permanent  thermocline  is  identified  at  100- 
250  m  depths  south  of  the  SPF  with  a  vertical  temper¬ 
ature  gradient  weakening  from  December  (near  0.1°C 
m  1 )  to  March  (near  0.05°C  m  1 ) .  The  permanent  ther¬ 
mocline  is  identified  at  25-100-m  depths  north  of  the 
SPF  in  December  with  a  vertical  gradient  near  0.06°C 
m  1 ,  which  is  much  weaker  than  that  south  of  the  SPF. 
From  January  to  March,  there  is  almost  no  evident  ther¬ 
mocline  north  of  the  SPF. 

During  the  prevailing  summer  monsoon  season  (Jun- 
Aug),  a  shallow  seasonal  thermocline  occurs  in  the 
whole  JES  basin  with  a  much  stronger  strength  north 
of  rather  than  south  of  the  SPF;  and  overlays  relatively 
uniform  water  north  of  the  SPF  and  stratified  water  (the 
permanent  thermocline)  south  of  the  SPF.  North  of  the 
SPF  a  seasonal  thermocline  appears  near  the  surface 
(above  50-m  depth)  with  a  vertical  gradient  enhancing 
from  0.25°C  m  1  in  June  to  0.36°C  m_I  in  August.  This 
strong  and  shallow  thermocline  isolates  the  exchange 
of  seawater  below  the  thermocline  from  atmospheric 
forcing  and  makes  this  water  (north  of  the  SPF  under 
the  thermocline)  quite  uniform.  South  of  the  SPF  a  sea¬ 
sonal  thermocline  is  wider  (25-100-m  depths)  and 
weaker  with  a  vertical  gradient  around  0. 13°C  m  1 .  Such 
a  north-south  asymmetric  pattern  in  July  was  previously 
presented  by  Kim  and  Kim  (1999)  using  the  Circulation 
Research  of  the  East  Asian  Marginal  Seas  (CREAMS) 
data  taken  mainly  in  the  summer  of  1995.  However,  the 
seasonal  thermocline  north  of  the  SPF  is  much  stronger 
from  the  GDEM  data  than  from  the  CREAMS -1995 
data. 

The  seasonal  thermocline  has  a  strong  asymmetry 
across  the  SPF.  North  of  the  SPF,  the  seasonal  ther¬ 
mocline  occurs  near  the  surface  in  April  and  May,  en¬ 


hances  drastically  during  the  summer  monsoon  season, 
survives  in  the  fall,  and  is  still  quite  strong  with  a  ver¬ 
tical  gradient  of  0.1 2°C  m  1  in  November.  It  weakens 
drastically  in  December.  South  of  the  SPF,  the  seasonal 
thermocline  occurs  in  the  summer  monsoon  season,  sur¬ 
vives  in  the  early  fall,  and  disappears  in  November. 

6.  Salinity 

a.  Sea  surface 

The  monthly  sea  surface  salinity  (SSS)  held  (Fig.  10) 
shows  a  strong  seasonal  variation  with  less  (more)  hor¬ 
izontal  variability  in  winter  (summer).  The  saline  Ku- 
roshio  water  (>34.2  psu)  enters  the  JES  from  the  Tsush¬ 
ima/Korean  Strait  into  the  JES  and  forms  two  permanent 
salty  centers  located  in  the  northern  JB  (west  of  Hok¬ 
kaido  Island)  with  the  salinity  higher  than  34.0  psu,  and 
the  area  between  UTB  and  YB  with  the  maximum  sa¬ 
linity  of  34.3  psu  in  June,  respectively.  The  northern  JB 
salty  center  has  less  seasonal  variation  than  the  UTB/ 
YB  salty  center.  Around  the  UTB/YB  salty  center,  there 
are  several  fresh  centers.  The  winter  (Feb)  field  (Fig. 
10)  is  consistent  with  that  reported  by  Kim  and  Kim 
(1999,  Fig.  9  in  their  paper)  using  the  dataset  of  the 
Japan  Oceanographic  Data  Center  during  1930-90. 

Three  important  features  are  found  from  the  monthly 
SSS  anomaly  held  (Fig.  11):  (i)  a  basinwide  positive 
(negative)  anomaly  occurs  in  the  winter  (summer);  (ii) 
the  seasonal  variability  of  the  surface  salinity  held  is 
stronger  in  the  south  than  in  the  north  of  the  SPF;  and 
(iii)  there  are  two  saline  activity  centers  located  at  Tartar 
Strait  and  Tsushima  Strait,  respectively. 

The  SSS  seasonal  variation  is  characterized  by  the 
transition  of  the  two  patterns:  the  winter  pattern  (bas¬ 
inwide  positive  anomaly  with  the  maximum  value  of 
0.8  psu  near  Tsushima  Strait  in  March)  and  the  summer 
pattern  (basinwide  negative  anomaly  with  the  minimum 
value  of  —0.8  psu  near  Tsushima  Strait  in  August-Sep- 
tember).  The  winter  pattern  prevails  in  January-March, 
and  the  summer  pattern  prevails  in  July-September.  The 
basinwide  positive  (negative)  SSS  anomaly  in  the  winter 
(summer)  might  be  caused  by  more  (less)  evaporation 
due  to  strong  winter  (weak  summer)  monsoon  winds 
and  less  (more)  precipitation. 

In  winter  (Jan-Mar),  two  salty  centers  occur  at  both 
Tartar  and  Tsushima  Straits,  respectively.  The  positive 
Tartar  Strait  SSS  anomaly  is  greater  than  0.2  psu  with 
a  maximum  value  of  0.4  psu  in  January  and  February. 
The  positive  Tsushima  Strait  SSS  anomaly  is  greater 
than  0.4  psu  with  a  maximum  value  of  0.8  psu  in  March. 
The  two  saltier  centers  are  caused  by  different  physical 
processes:  brine  ejection  due  to  ice  freezing  for  the  Tar¬ 
tar  saltier  center,  and  the  intrusion  of  open  ocean  water 
for  the  Tsushima  saltier  center. 

The  spring  (Apr-May)  is  characterized  by  the  oc¬ 
currence  and  southward  expansion  of  a  fresher  (negative 
anomaly)  center  in  Tartar  Strait  with  a  minimum  value 
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Fig.  8.  Zonal  cross  section  of  the  monthly  mean  temperature  (°C)  at  (a)  37°  and  (b)  43°N.  Contour  interval:  1°C. 
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Fig.  8.  (Continued) 
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Fig.  9.  Latitudinal  cross  section  (135°E)  of  the  monthly  mean  temperature  (°C).  Contour  interval:  1°C. 
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Fig.  10.  Monthly  mean  salinity  (psu)  field  at  the  ocean  surface.  Low  salinity  values  (5  <  34.06  psu)  are 

shaded. 
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Fig.  11.  Monthly  mean  salinity  anomaly  (psu)  field  at  the  ocean  surface.  Negative  values  are  shaded. 
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of  —0.4  psu.  Freshwater  influx  due  to  ice  melting  in  the 
spring  is  responsible  for  the  Tartar  Strait  negative  anom¬ 
aly. 

In  June,  the  negative  salinity  anomaly  occupies  north 
of  the  SPF,  and  another  fresher  center  is  formed  near 
Tsushima  Strait  with  a  minimum  value  of  —0.6  psu 
south  of  Tsushima  Strait.  The  physical  process  causing 
this  negative  SSS  anomaly  is  that  in  summer  the  Tsush¬ 
ima  Current  carries  in  a  large  amount  of  Yellow  Sea 
water  that  is  diluted  by  river  runoff  during  the  summer 
monsoon  (Tomczak  and  Godfrey  1994).  The  northward 
expansion  of  this  fresher  anomaly  into  the  JES  basin 
leads  to  the  summer  SSS  pattern. 

In  summer  (Jul-Sep),  a  strong  fresher  center  with  a 
minimum  value  of  —1.0  psu  near  Tsushima  Strait  ex¬ 
pands  northward  and  sweeps  the  whole  JES  with  neg¬ 
ative  anomalies.  The  SSS  anomaly  is  around  —0.2  psu 
north  of  the  SPF  and  from  —0.2  to  —1.0  south  of  the 
SPF. 

The  autumn  (Oct-Nov)  is  characterized  by  the  oc¬ 
currence  of  two  positive  anomaly  areas:  one  near  Tartar 
Strait  and  the  other  near  Tsushima  Strait.  The  expansion 
of  the  two  positive  anomaly  areas  into  the  JES  basin 
leads  to  the  transition  to  the  winter  pattern. 

b.  Intermediate  level  (150  m) 

Figure  12  shows  horizontal  depictions  of  the  monthly 
mean  salinity  at  150-m  depth  with  a  contour  interval  of 
0.1  psu.  The  depth  of  150  m  corresponds  to  the  middle 
level  of  the  JIW  (low  salinity)  as  well  as  the  HSIW 
(high  salinity).  The  low  salinity  water  of  the  JIW  less 
than  34.06  psu  (shaded  areas  in  Fig.  12)  prevails  in  the 
western  JB  and  the  UTB.  It  expands  eastward  during 
the  summer  monsoon  season  (Jun-Sep)  and  retreats 
westward  during  the  winter  monsoon  season  (Dec- 
Mar).  The  high  salinity  water  of  the  HSIW  occupies  the 
eastern  JB.  This  is  consistent  with  earlier  work  by  Kim 
and  Kim  (1999).  The  saline  water  occupies  the  east  and 
south  and  less  saline  water  appears  in  the  north  and 
west,  especially  at  the  northwestern  boundary.  The  34. 1 
isohaline  is  collocated  with  the  SPF  (Fig.  6)  with  the 
salinity  above  (below)  34.1  in  the  south  (north)  of  the 
SPF.  North  of  the  SPF  salinity  is  relatively  uniform 
(34.0-34.1  psu)  throughout  the  year.  South  of  the  SPF 
salinity  changes  from  34.1  to  34.5  psu. 

Figure  13  shows  horizontal  depictions  of  the  monthly 
mean  salinity  anomaly  at  150-m  depth  with  a  contour 
interval  of  0.05  psu.  The  absolute  values  of  the  anomaly 
are  less  than  0.2  psu  throughout  the  JES  basin.  The 
winter  (Jan-Mar)  pattern  is  featured  by  a  major  positive 
anomaly  north  of  the  SPF  and  negative  anomaly  south 
of  the  SPF.  The  summer  (Jul-Sep)  pattern  seems  a  mirror 
image  of  the  winter  pattern:  a  major  negative  anomaly 
north  of  the  SPF  and  positive  anomaly  south  of  the  SPF. 
The  spring  and  fall  anomalies  are  transitions  between 
the  two  patterns. 


c.  Zonal  cross  sections  (37°  and  43  °N) 

Zonal  cross  sections  (37°  and  43°N)  of  the  monthly 
mean  salinity  show  an  evident  SMIN  south  of  the  SPF 
(Fig.  14a)  and  absence  of  a  SMIN  north  of  the  SPF  (Fig. 
14b).  This  is  consistent  with  many  earlier  studies  such 
as  Miyazaki  (1952,  1953),  Miyazaki  and  Abe  (1960), 
Kim  and  Chung  (1984),  Senjyu  (1999),  and  Kim  and 
Kim  (1999). 

South  of  the  SPF  at  37°N  (Fig.  14a),  a  strong  (sea¬ 
sonal)  halocline  occurs  from  the  surface  to  depth  30  m 
in  June  (0.01  psu  m~‘),  intensifies  during  the  summer 
monsoon  season  to  a  maximum  value  of  around  0.03 
psu  m  1  in  August,  and  weakens  from  September  to 
December.  In  January,  the  halocline  erodes  and  disap¬ 
pears.  During  the  summer  monsoon  season  (Jun-Sep), 
the  SMIN  shifts  eastward  from  the  UTB  (131°-133°E) 
in  June  to  the  southern  YB  (1320-135°E)  in  July-Sep- 
tember.  A  horizontally  oriented  SMAX  (.S'  >  34.2  psu) 
appears  above  the  SMIN  with  the  interface  at  depths 
200-300  m.  In  October,  this  horizontally  oriented 
SMAX  is  broken  into  two  parts  with  each  part  enclosed 
by  the  34.2  psu  isohaline.  A  narrow  relatively  low  sa¬ 
linity  area,  sandwiched  by  the  two  SMAX  centers,  oc¬ 
curs  in  the  central  UTB  (near  132°E)  beneath  depth  100 
m  and  strengthens  from  November  to  January  accom¬ 
panying  the  erosion  of  the  surface  halocline.  In  January, 
the  near-surface  halocline  disappears,  the  SMIN  (en¬ 
closed  by  34.1  psu  isoline)  appears  in  the  UTB  (131°- 
133°E)  and  reaches  the  surface,  and  the  two  SMAX 
centers  weaken  (especially  the  western  one).  The  SMIN 
expands  horizontally,  descending  vertically,  and  be¬ 
comes  very  evident  in  February.  From  February  to  June, 
the  SMIN  continues  to  descend.  The  top  of  the  SMIN 
occurs  at  depths  150-200  m  in  February  and  at  200- 
300  m  in  April-June. 

North  of  the  SPF  at  43°N  (Fig.  14b),  the  SMIN  shows 
up  in  the  upper  layer  (above  100  m)  of  the  western  JB 
(west  of  136°E)  year-round.  This  is  consistent  with  Kim 
and  Kim’s  (1999)  identification  that  the  JIW  in  the  west¬ 
ern  JB  is  characterized  by  low  salinity  ( S  <  34.06  psu). 
The  high  salinity  water  ( S  >  34.1  psu)  of  the  HSIW 
appears  in  the  central  (136°-138°E)  and  the  eastern 
(138.5o-140°E)  JB.  The  eastern  JB  SMAX  appears  al¬ 
most  year-round  except  during  the  spring  season  (Mar- 
May).  The  central  JB  SMAX  is  not  as  evident  as  the 
eastern  one  except  in  February  and  March. 

d.  Latitudinal  cross  section  (135°E) 

The  strong  north-south  haline  asymmetry  across  the 
SPF  is  also  identified  from  a  latitudinal  cross  section 
(135°E)  of  the  monthly  mean  salinity:  the  appearance 
(disappearance)  of  the  SMIN  south  (north)  of  the  SPF 
(Fig.  15). 

South  of  the  SPF,  SMIN  ( S  <  34.06  psu)  occurs  dur¬ 
ing  the  summer  monsoon  season  (Jul-Oct)  underneath 
a  horizontally  oriented  SMAX  with  a  salty  core  ( S  > 
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Fig.  12.  Monthly  mean  salinity  (psu)  field  at  150-m  depth.  Low  salinity  values  ( S  <  34.06  psu)  are  shaded. 
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Fig.  13.  Monthly  mean  salinity  anomaly  (psu)  field  at  150-m  depth.  Negative  values  are  shaded. 
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Fig.  14.  Zonal  cross  section  of  the  monthly  mean  salinity  (psu)  at  (a)  37°  and  (b)  43°N.  Contour  interval:  0.1  psu. 
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Fig.  14.  ( Continued ) 
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Fig.  15.  Latitudinal  cross  section  (135°E)  of  the  monthly  mean  salinity.  Contour  interval:  0.1  psu. 
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34.5  psu)  at  depth  100  m.  The  SMAX  with  a  salty  core 
(, S  =  34.3  psu)  was  observed  in  October  1969  when  a 
hydrographic  survey  for  the  whole  JES  basin  was  car¬ 
ried  out  by  the  Japan  Meteorological  Agency  (Kim  and 
Kim  1999).  The  interface  between  the  SMIN  and  the 
SMAX  is  located  at  depth  200-300  m.  During  the  winter 
monsoon  season,  the  SMIN  is  not  evident. 

North  of  the  SPF,  high-salinity  water  ( S  >  34. 1  psu) 
of  the  HSIW  appears  in  the  central  JB  (40°-43°N)  with 
strengthening  in  the  winter  and  weakening  in  the  sum¬ 
mer. 

7.  Absolute  geostrophic  velocity 

The  inverted  absolute  geostrophic  velocity  field  (Fig. 
16)  coincides  with  earlier  descriptions  of  the  JES  cir¬ 
culation  (e.g.,  Uda  1934;  Suda  and  Hidaka  1932;  Suda 
et  al.  1932;  Suk  1988;  Kim  and  Yoon  1996).  However, 
the  currents  at  the  four  major  straits  are  not  well  in¬ 
verted.  For  example,  the  velocity  at  the  Tsushima/Korea 
Strait  is  too  small,  and  there  is  no  velocity  inverted  at 
the  other  straits  such  as  Tartar  Strait.  This  may  be  caused 
by  the  narrowness  and  shallowness  of  the  straits. 

The  inverted  absolute  geostrophic  velocity  field 
shows  an  evident  branching  of  the  TWC.  North  of  35°N, 
the  TWC  bifurcates  into  an  eastern  channel  (first  branch) 
and  a  western  channel.  The  strength  of  the  TWC  at  both 
channels  reduces  with  depth.  Flow  through  the  western 
channel  (i.e.,  EKWC)  closely  follows  the  Korean  coast 
until  it  separates  near  37°-38°N  into  two  branches:  the 
offshore  branch  (or  the  second  branch)  follows  the  Polar 
Front  to  the  western  coast  of  Hokkaido  Island  and  the 
alongshore  branch  (or  the  third  branch)  moves  north¬ 
ward  (i.e.,  the  EKWC).  Such  a  three-branch  pattern  was 
first  identified  by  Suda  and  Hidaka  (1932)  and  Suda  et 
al.  (1932)  using  hydrographic  and  current  meter  data 
from  June  to  September  in  1929  and  in  1930,  respec¬ 
tively.  Since  then,  the  existence  of  three  branches  has 
been  believed  to  be  the  typical  TWC  flow  pattern.  Using 
temperature  and  salinity  data  mainly  obtained  in  1973, 
Kawabe  (1982a)  investigated  the  seasonal  variation  of 
the  TWC  branching  and  found  that  the  first  branch  exists 
in  spring  and  summer,  the  second  branch  only  in  summer 
from  June  to  August,  and  the  third  branch  in  all  seasons. 
His  observational  results  were  simulated  using  a  two- 
layer  numerical  model  (Kawabe  1982b;  Sekine  1986). 

The  inverted  flow  through  the  eastern  channel  (first 
branch),  which  closely  follows  the  Japanese  coast  (i.e., 
JNB),  is  weaker  than  the  flow  through  the  western  chan¬ 
nel  (Fig.  16).  The  maximum  speed  of  the  first  branch 
is  found  around  0.1  m  s  1 ,  occurring  near  the  Japanese 
coast  (35.5°N,  135°E)  at  the  surface.  The  JNB  is  too 
weak  to  be  identified  in  the  winter  (Dec-Feb).  In  March, 
the  JNB  starts  along  the  Japanese  coast  and  strengthens 
in  the  spring  and  reaches  its  maximum  velocity  in  July. 
This  coincides  with  Kawabe’s  (1982a)  observational  re¬ 
sults. 

The  inverted  flow  along  the  Polar  Front  (second 


branch)  may  be  called  the  Polar  Front  Current  (PFC). 
The  maximum  speed  of  the  PFC  is  found  at  the  surface 
around  0.2  m  s_1  (Fig.  16).  An  asymmetric  cyclonic 
gyre  is  found  in  the  mid-JES  (38°-44°N,  135°-140°E) 
in  the  winter  monsoon  season  (Dec-Mar)  with  the  flow 
associated  with  the  PFC  at  the  southern  and  eastern 
flanks  and  the  flow  from  the  north  (mid-JES  cold  cur¬ 
rent)  at  the  northern  and  western  flanks.  The  seasonal 
variability  of  this  gyre  is  largely  determined  by  the  sea¬ 
sonal  variability  of  the  PFC.  The  PFC  has  a  weak  sea¬ 
sonal  variation  in  flow  pattern  and  a  strong  seasonal 
variation  in  current  speed. 

An  interesting  feature  is  the  out-of-phase  variation 
between  the  PFC  (along  the  Polar  Front)  and  the  JNB 
(along  the  west  coast  of  Japan).  For  example,  the  PFC 
(JNB)  strengthens  (weakens)  from  July  to  September, 
and  the  PFC  (JNB)  weakens  (strengthens)  in  spring  from 
March  to  April. 

The  seasonal  variation  of  the  JES  major  currents  is 
also  shown  in  the  latitudinal  cross  section  of  the  monthly 
mean  u  component  (Fig.  17).  Two  eastward  flowing 
currents,  the  JNB  (near  36°N)  and  the  PFC  (38°-40°N), 
are  well  represented.  The  PFC  has  the  maximum  speed 
of  0.22  m  s_1  (minimum  speed  of  0.12  m  s_l)  occurring 
in  September  (Feb-Apr)  at  the  surface.  During  the  sum¬ 
mer  monsoon  season  (Jul-Sep),  the  PFC  strengthens 
(maximum  speed  increases  from  0.18  m  s_1  to  0.22  m 
s  1 ),  and  the  JNB  weakens  (maximum  speed  decreases 
from  0.12  m  s_I  to  0.04  m  s_1).  During  the  winter  mon¬ 
soon  season  (Dec-Apr),  the  PFC  weakens  (maximum 
speed  decreases  from  0.18  m  s_1  to  0.12  m  s_1),  and 
the  TWC  strengthens  (maximum  speed  increases  from 
0.02  m  s_l  in  January  to  0.1  ms1  in  April).  The  PFC 
volume  transport  shows  an  evident  seasonal  variation 
with  the  maximum  value  of  1.32  Sv  (Sv  =  106  m3  s_1) 
in  November-December  and  the  minimum  value  of  0.88 
Sv  in  April  (Fig.  18). 

8.  Conclusions 

The  goal  of  this  study  is  to  investigate  the  seasonal 
variabilities  of  the  thermohaline  structure  and  the  cir¬ 
culation  patterns  of  the  JES  using  the  public  domain 
GDEM  climatological  hydrographic  data.  The  major  re¬ 
sults  from  this  research  are  summarized  as  follows: 

1)  The  seasonal  thermal  variability  is  evident  at  the 
surface  and  weakens  with  depth.  However,  the  SPF  ex¬ 
ists  continuously  through  the  year.  Differential  warm¬ 
ing/cooling  is  found  across  the  SPF.  At  the  intermediate 
level  (depth  150  m),  the  SPF  is  situated  almost  at  the 
same  location  as  at  the  surface.  The  winter  cooling 
(summer  warming)  is  much  stronger  south  of  the  SPF 
(anomaly  larger  than  2°C)  than  north  of  the  SPF  (anom¬ 
aly  less  than  1°C).  Such  a  differential  cooling  (warming) 
causes  the  northward  increase  (decrease)  of  temperature 
anomaly  across  the  SPF,  which  implies  the  reduction 
(enhancement)  of  the  SPF  strength  in  the  winter  (sum¬ 
mer). 
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Fig.  16.  Inverted  monthly  mean  surface  circulation. 
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Fig.  17.  Latitudinal  cross  sections  of  the  inverted  monthly  mean  u  velocity  along  135°E.  Contour  interval:  2  cm  s 
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Fig.  18.  Monthly  variation  of  the  volume  transport  (Sv)  of  the 
Polar  Front  Current  at  135°E  cross  section. 

2)  The  seasonal  thermal  variability  has  north-south 
asymmetry  with  a  strong  seasonal/weak  permanent  ther- 
mocline  structure  south  of  the  SPF  and  a  strong  sea¬ 
sonal/weak  permanent  thermocline  structure  north  of  the 
SPE  South  of  the  SPF,  the  permanent  thermocline  is 
located  at  depths  80-125  m  appearing  year-round  with 
the  maximum  strength  (0.12°C  m  1 )  in  August.  Above 
it,  the  seasonal  thermocline  occurs  from  the  surface  to 
depth  50  m  in  June  (0.15°C  m-1),  intensifies  during  the 
summer  monsoon  season  to  a  maximum  value  of  around 
0.36°C  m  1  in  August,  and  weakens  in  September.  In 
October,  the  seasonal  thermocline  erodes  and  the  ocean 
mixed  layer  forms  and  then  deepens  in  winter.  The  ocean 
mixed  layer  starts  to  warm  at  a  rate  of  2°C  month-1 
from  March  to  May,  and  its  depth  shoals. 

North  of  the  SPF,  the  permanent  thermocline  is  quite 
weak  and  located  near  the  surface  to  depth  300  m  with 
an  evident  westward  uplift  from  the  Japanese  coast 
(thickness  around  300  m)  to  the  Russian  coast  (thickness 
around  100  m).  The  seasonal  thermocline  occurs  from 
the  surface  to  depth  50  m  in  May  (~0.08°C  m-1),  in¬ 
tensifies  during  the  summer  monsoon  season  to  a  max¬ 
imum  value  of  around  0.5°C  m  1  in  August  and  Sep¬ 
tember,  and  weakens  in  October.  In  November,  the  sea¬ 
sonal  thermocline  erodes  and  becomes  the  part  of  the 
permanent  thermocline,  which  weakens  during  the  pre¬ 
vailing  winter  monsoon  season. 

3)  Three  important  features  are  found  from  the 
monthly  sea  surface  salinity  anomaly  field:  a  basinwide 
positive  (negative)  anomaly  occurs  in  winter  (summer); 
the  seasonal  variability  of  the  surface  salinity  field  is 
stronger  south  than  north  of  the  SPF;  and  there  are  two 
saline  activity  centers  located  at  Tartar  Strait  and  Tsush¬ 
ima  Strait,  respectively. 

4)  The  seasonal  haline  variability  has  a  strong  north- 
south  asymmetry  with  an  evident  salinity  minimum 
south  of  the  SPF  and  absence  of  a  salinity  minimum 
north  of  the  SPF.  South  of  the  SPF,  the  SMIN  ( S  < 
34.06  psu)  occurs  during  the  summer  monsoon  season 
(Jul-Oct)  underneath  a  horizontally  oriented  SMAX 
with  a  salty  core  ( S  >  34.5  psu)  at  depth  100  m.  The 


interface  between  the  SMIN  and  the  SMAX  is  located 
at  depth  200-300  m.  During  the  winter  monsoon  season, 
the  SMIN  becomes  less  evident.  North  of  the  SPF,  the 
high,  salinity  water  ( S  >  34.1  psu)  of  the  HSIW  appears 
in  the  central  JB  (40°-43°N),  strengthening  in  winter 
and  weakening  in  summer. 

5)  The  inverted  absolute  geostrophic  velocities  rep¬ 
resent  the  JES  circulation  reasonably  well,  especially 
the  Tsushima  Current  and  its  bifurcation.  The  Tsushima 
Current  bifurcates  into  a  western  and  an  eastern  branch 
north  of  the  Korea/Tsushima  Strait.  The  maximum  speed 
of  the  western  branch  is  found  around  0.2  ms-1,  oc¬ 
curring  near  the  SPF  at  the  surface.  The  maximum  speed 
of  the  eastern  branch  is  found  around  0.1  m  s_1,  oc¬ 
curring  near  the  Japanese  coast  (135°E,  35.5°N)  at  the 
surface.  The  strength  of  the  Tsushima  Current  in  both 
branches  reduces  with  depth.  An  evident  cyclonic  gyre 
is  found  in  the  mid-JES  (38°-44°N,  135°-140°E)  with 
the  flow  associated  with  the  SPF  at  its  eastern  flank  and 
the  flow  from  the  north  (mid-JES  cold  current)  at  its 
western  flank. 

6)  The  SPF  current  has  a  weak  seasonal  variation  in 
flow  pattern,  a  strong  seasonal  variation  in  current 
speed,  and  out-of-phase  variability  with  the  Tsushima 
Warm  Current.  When  the  SPF  current  strengthens  from 
July  (0.18  m  s-1)  to  September  (0.22  m  s-1),  the  Tsush¬ 
ima  Warm  Current  weakens  from  0.12  m  s  1  to  0.04  m 
s_1.  When  the  SPF  current  weakens  from  January  (0.18 
m  s-1)  to  April  (0.12  m  s-1),  the  Tsushima  Warm  Cur¬ 
rent  strengthens  from  0.02  m  s-1  to  0.1  m  s  1 .  The 
seasonal  variability  of  the  mid-JES  cyclonic  gyre  is 
largely  determined  by  the  seasonal  variability  of  the  SPF 
current.  The  volume  transport  of  the  SPF  current  has  a 
maximum  value  of  1.32  Sv  in  November-December  and 
a  minimum  value  of  0.88  Sv  in  April. 
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